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Techniques for detesnnlning glucose levels 



5 Technical Field 

The invention relates devices for the deter- 
mination of glucose, to methods for calibrating or oper- 
ating such devices and to methods for measuring glucose. 

10 

Background Art 

It has been known that the glucose of living 
15 tissue can be measured non-invasively by applying a sen- 
sor arrangement, in particular an electrode arrangement 
to the skin of a patient and measuring the response of 
the electrode arrangement to a suitable electric signal. 
Such a technique is described in WO 02/069791, the dis- 
20 closure of which is enclosed herein in its entirety. 

Even though this type of device is well able 
to monitor glucose, it needs careful calibration accord- 
ing to a well-defined protocol and must be operated under 
defined conditions in order to yield results of high ac- 
25 curacy. 

One important purpose of these devices is to 
provide a prediction of the time when a patient's glucose 
level may exceed certain limits. In particular, an early 
prediction of a possible hypoglycemia or hyperglycemia is 
30 desirable such that a patient or accompanying person may 
take adequate steps to prevent such a state. 



Disclosure of the Invention 

35 

Hence, in a first aspect, it is an object of 
the invention to provide a device and method of the type 
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mentioned above that allow a more accurate measurement of 
the glucose level in a living body. 

This object is achieved by the device and 
method of claims 1 and 26. 
5 According to this aspect of the invention, at 

least two temperatures are measured, wherein the first 
temperature depends in different manner on the skin tem- 
perature of the body and on the environmental temperature 
than the second temperature when the device is mounted on 

10 the body in its position of operation. This allows to 
compensate for the influence of both, the skin tempera- 
ture and the environmental temperature, which has been 
found to be advantagous because, as discussed in the de- 
tailed description, both temperatures affect the signals 

15 measured with the sensor arrangement in different man- 
ners . 

In a second aspect, the object of the first 
aspect is achieved by improving the calibration of the 
device. Two different calibration mechanisms are sug- 
20 gested, which can be used alternatively or in combina- 
tion. 

In both these mechanisms, the device is as- 
sumed to calculate, in normal operation, the glucose 
level from a function of the type F(si, .S2, ... Sjjj, slq, 

25 ai, ... ajj) , where F depends on input values si ... sj^ 
and calibration parameters ag ... ajy^. 

In a calibration phase, series of reference 
values g(ti) are obtained at times t^, e.g. using conven- 
tional glucose measurements. In the same phase, a series 

30 of raw input values sj(t'i) are measured at times t'^, 
which generally do not necessarily coincide with the 
times ti . At least part of the parameters slq ... a-^ are 
then derived from these measurements by coitvparing the 
values obtained by function F against the reference val- 

35 ues or against values derived from the reference values. 

In most cases, the number of times the input 
values s(t'i) have been measured will be considerably' 
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larger than the number of reference values gCt^) . Hence, 
in the first mechanism, in order to fully exploit all 
data, a prediction (interpolation) of the glucose level g 
at the times t'i is calculated from the reference values 
5 g{ti) Then the deviation of the values calculated by 

function F for the input values Sj(t'i) and the predicted 
glucose levels at times t' j is minimized by varying the 
parameters ag - . . ajj^/ thereby obtaining a set of cali- 
brated parameters . 

10 In the second mechanism, a "^^shift correction" 

is carried out during the calibration phase. For this 
purpose, the times are detected at which the device 
has shifted or moved in relation to the body during cali- 
bration. Such shifts generally cause the measured signals 

15 to change. When comparing the values obtained by fiinction 

F against the reference values or against values derived 

from the reference values as mentioned above, at least 

. one parameter ag is replaced by a sum 
P 

S aoi • jbi(t) 
i=0 

20 

with bi(t) being 1 (or, equivalently , any other non-zero 
constant value) f or .T^ < t < Ti+i and 0 otheirwise. As ex- 
plained in the' detailed description, this allows to com- 
pensate for the effects of the shifts. 

25 In any case it may be advantageous to carry 

out a recalibration step, e.g. each time after putting on 
the device. In this step, one of the parameters is re- 
calibrated to find axi optimum agreement between the glu- 
cose level calculated from the function F and a glucose 

30 level from a reference measurement. The reference meas- 
urement can e.g. be a conventional measurement, such as 
an invasive measurement. This allows to compensate for an 
offset caused by removing and remoimting the device. 

In another aspect, an object of the invention 

35 is to provide a device and method that are capable to 
provide early and reliable prediction of a possible hy- 



wo 2005/053526 



4 



PCT/IB2003/005704 



per-, or hypoglycemia. This object is achieved by the de- 
vice and method of claims 16 and 37 . 

The various aspects and mechanisms can be 
used in combination or separately. 

5 

Brief Description of the Drawings 

The invention will be better understood and 
10 objects other than those set forth above will become ap- 
parent when consideration is given to the following de- 
tailed description thereof . Such description makes refer- 
ence to the annexed drawings, wherein: 

Fig. 1 is a cross section of a device for 
15 measuring a glucose level. 

Fig. 2 is a block circuit diagram of the de- 
vice of Fig. 1, 

Fig. 3 is an apparatus for calibrating the 

device, 

20 Fig. 4 illustrates an advantageous aspect of 

the calibration method. 

Fig. 5 shows a signal shift upon a displace- 
ment of the device, and 

Fig. 6 illustrates a worst-case prediction of 
25 glucose levels with and without limits for the second or- 
der derivative. 



Modes for Carrying Out the Inv e ntion 

30 

The device: 

Fig. 1 shows a cross section of a device 100 ' 
for measuring a patient's glucose level. It comprises a 
housing 1 closed on one side by an electrode plate 2 . A 
35 display 3 is arranged opposite electrode plate 2 . Elec- 
tronic circuitry is arranged between electrode plate 2 
and display 3 . 
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Electrode plate 2 comprises an electrically 
insulating substrate 4 . A strip electrode 5 covered by an 
insulating layer 5a and a ring electrode 6 are arranged 
on an outer side 7 of insulating substrate 4 . An inner 
5 side 8 of insulating substrate 4 is covered by a ground 
electrode 9. A plurality of though- contacts 10 connect 
ring electrode 6 to ground electrode 9. A further 
through-contact 11 connects one end of strip electrode 5 
to a contact pad 12 arranged on inner side 8. 

10 A first teitiperature sensor 15 is mounted to 

gro\ind electrode 9 in direct thermal contact thereto . The 
large number of through-contacts 10 ensures that groiind 
electrode 9 follows the temperature of ring electrode 6 
and therefore the temperature of the specimen, the sur- 

15 face of which is indicated by a dotted line 16, closely. 

Leads 18 are provided to connect ground elec- 
trode 9, contact pad 12 and first temperature sensor 15 
to the electronic circuitry arranged on a printed circuit 
board 19 f oirming an assembly of electronic components . 

20 Printed circuit board 19 is advantageously arranged on a 
side of the device that is substantially opposite to the 
side of electrode plate 2. A battery 21 for powering the 
circuitry is arranged between printed circuit board 19 
and electrode plate 2 . 

25 A second temperature sensor 22 is arranged on 

printed circuit board 19 and in direct thermal contact 
thereto . 

The design of the electrodes 5, 6, 9 of the 
present sensor can correspond to the one described in 
30 reference to Figs. 2 and 4 of WO 02/069791, which de- 
scription is enclosed by reference herein. 

Fig. 2 shows a block circuit diagram of the 
circuitry of device 100. It comprises a voltage con- 
trolled oscillator (VCO) 31 as a signal source for gener- 
35 ating a sine wave signal or another periodic signal. This 
signal is fed to two amplifiers 32, 33. The output of 
first amplifier 32 is connected via a resistor Rl to a 
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first signal path 34. A resonant circuit 35 comprising an 
inductance L and a capacitor C in series is connected be- 
tween first signal path 34 and ground. The output of sec- 
ond amplifier 33 is connected via a resistor R2 to a sec- 
5 ond signal path 3-6. Second signal path 36 can be substan- 
tially identical to first signal path 34 but contprises a 
resistor R3 as a reference load instead of resonant cir- 
cuit 35. 

Both signal paths 34, 36 are fed to a measur- 

10 ing circuit 37, which deteinnines the relative amplitude A 
of both signals and/or their mutual phase shift phi. 
Relative amplitude A can e.g. be the amplitude of first 
signal path 34 in units of the ait^litude of second signal 
path 3 6 (wherein the amplitudes are the peak values of 

15 the sine waves) . 

The output signal of measuring circuit 37 is 
fed to a microprocessor 38, which also controls the op- 
eration of VCO 31- 

Microprocessor 38 further samples the first 

20 and second temperature signals Tl, T2 from first and sec- 
ond temperature sensors 15, 22. It also controls display 
device 3, an input device 40 with user operable controls, 
and an interface 41 to an external computer. A memory 42 
is provided for storing calibration parameters, measure- 

25 ment results, further data as well as firmware for micro- 
processor 38. At least part of memory 42 is non-volatile. 

Inductance L of the device of Fig. 2 can be 
generated by a coil and/ or by the leads and electrodes of 
capacitor C- Its value is generally known with reasonable 

30 accuracy. 

Capacitor C of the device of Fig. 2 is foarmed 
between strip electrode 5 and ring electrode 6 and is 
used for probing the specimen. For this puirpose, the 
electrodes are arranged on the skin 16 of the patient as 
35 shown in Fig. 1. 

For a good and permanent contact with the pa- 
tient 's skin, the device is advantageously worn on an arm 
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or leg and provided with a suitable holder or wrist band 
43. 

The geometry of the electrodes is selected 
such that the electric field generated by them extends 
5 into the specimen and the body liquid to be measured. Ad- 
vantageously, at least one of the electrodes of the ca- 
pacitor is electrically insulated such that capacitor C 
is primarily a capacitive load, the capacitance and loss 
of which depend on the electrical properties (i.e. the 
10 response) of the specimen at the frequency of VCO 1. 

In summary, the device shown in Figs. 1 and 2 

comprises : 

- an electrode arrangement or sensor arrange- 
ment comprising the electrodes 5 and 6 and 

15 - processing circuitry including the elements 

31-33, 37, 38 for measuring the response of the sensor 
arrangement or electrode arrangement to an electrical 
signal and deriving the glucose level therefrom. 

In addition, it can comprise at least two 

20 temperature sensors 15, 22, the signals of which depend 
in different manner on the skin temperature of the body 
and on the environmental temperature. Both these teitpera- 
tures can be taken into account when determining the glu- 
cose level . 

25 * . 

Basic principle of operation: 
The basic principle of operation of the de- 
vice is described in WO 02/069791. 

To measure the concentration of glucose in 

30 the body fluid of the patient, microprocessor 38 can e.g. 
initiate a measurement cycle consisting of a frequency 
sweep of VCO 1. The sweep should start at a frequency 
fmax above the expected resonance frequency fO of the 
resonant circuit 5 and extend to a frequency fniin below 

35 resonance frequency f 0 (or vice versa) . Typical frequen- 
cies are given in WO 02/069791. During this sweep, the 
electrical properties of the two signal paths 34, 36 will 
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vary in different manner. The amplitude determined by 
measuring circuit A will fall to a minimum AO at a char- 
acteristic frequency fO, as described in WO 02/069791. At 
the- same time or close thereto, phase shift phi crosses 
5 zero. 

Microprocessor 38 measures AO and/ or fO as 
input values describing the physiological state of the 
patient's tissue- In addition to the input values of AO 
and/or fO, microprocessor 38 measures the temperature 

10 values Tl and T2 as further input values . Using suitable 
calibration data, the glucose level can be derived from 
these input values . 

Such calibration data can be determined in 
straightforward manner using methods Imown to the person 

15 skilled in the art. In the following, however, some ad- 
vantageous techniques are presented for the determination 
of the glucose level with the type of device described 
here as well as for its calibration. 

In general, microprocessor 38 will use a for- 

20 mula of the type 

g = F(si, S2, ... sjg, SLQ, ai, ... a^) (D 

for determining the glucose level g (or a parameter in- 
25 dicative thereof) from N measured input values s^, 32, 
... sjj (N > 0) , where the function F has M+1 parameters 
SQ/ ^If ••• (M > 0)/ at least some of which have to be 
determined in suitable calibration experiments. 

The measured input values Si are e.g. values 
30 directly or indirectly derived from the amplitude AO, the 
corresponding frequency fO, and the temperatures Tl, T2 . 
The input values can e.g. be the most recent values meas- 
ured or they can be a time average or a median over a 
given number of recent measurements. 
35 In an advantageous embodiment, the values s^ 

= AO, S2 = fO, S3 = Tl and S4 = T2 are used. 
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The fimction F can be empirical or it can be 
based at least partially on a model describing the physi- 
cal nature of the mechanisms involved. 

Under the approximation that the relation be- 
5 tween the glucose level g and the measured values Si is 
linear, we have 

g = ao + ai-si + a2-S2 + ... au-s^ (2a) 

with M = N. 

10 Equation (2a) has the advantage of being lin- 

ear in the input values si as well as the parameters aj , 
which simplifies calibration as well as evaluation. More 
refined models can, however, be used as well. 

15 Temperature compensation: 

It is presently understood that electrical* 
properties of the topmost skin layers and therefore of 
the signals AO and f 0 depend not only on the glucose 
level, but also on the temperature Ts of the skin and un- 

20 derlying tissue as well as on the temperature Te of the 
environment. This is at least in part due to the fact 
these properties depend on the amoimt of blood in the 
skin and underlying tissue, which in turn affects the 
temperature of the skin. Hence, it is advantageous to 

25 ' measure the skin temperature, a first approximation of 
which can be derived from the signal from temperature 
sensor Tl. 

However, the skin temperature is not only a 
function of the amount of blood in the skin and underly- 

30 ing tissue, but also of the environmental temperature Te. 
Hence, it is also advantageous to measure the environ- 
mental temperature, a first approximation of which can be 
derived from the signal from temperature sensor T2 . 

Hence, device 100 is advantageously equipped 

35 with at least two temperature sensors Tl and T2 , the sig- 
nals of which depend in different manner on the tempera- 
tures Ts and Te, such that a measurement of Tl and T2 is 
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indicative of both temperatures Ts and Te. Hence, at 
least one of the input values s^ should be derived from 
the signal of first temperature sensor 15 and at least 
another one of the input values s^ should be derived from 
5 the signals of second temperature sensor 22 • 

Advantageously, one of the temperature sen- 
sors is closer to the electrodes 5, 6 (and therefore to 
the body of the patient) than the other sensor. For exam- 
ple, the first temperature sensor 15 is arranged at the 

10 same side of housing 1 as the electrodes 5, 6 and the 
second temperature sensor 22 at the opposite side. 

The measured values may also depend on the 
temperature of the electronic circuits because the prop- 
erties of voltage sources, A/ D- converters and other cir- 

15 cuitry are generally temperature dependent. Hence, it may 
also be advantageous to measure a temperature that is in- 
dicative of the circuit temperature Tc. In the present 
embodiment, this is especially true for temperature T2, 
i.e. by using the signal from second teinperature sensor 

20 22, changes of the circuit temperature Tc can be ac- 
counted for. However, an additional third temperature 
sensor for specifically measuring circuit temperature Tc 
may be provided as well. 



25 Calibration: 

In the following, advantageous methods for 
calibrating the device are described. 

a) Parameter Determination: 
30 A basic calibration of the device is required 

for each new patient - 

In a first step of the basic calibration, the 
patient undergoes a calibration phase in which the glu- 
cose level is measured repetitively by an alternative 
35 method of measurement, e.g. by a conventional invasive 
technique, in order to obtain a series of K reference 
values g(ti), g(t2)/ gCt^) at times t^ through t^- In 
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the same period, the input values are measured repeti- 
tively at L times t'l through t'L, wherein L can be much 
larger than K. All measured values Si(t'j) (i = 1 ... N, 
j =1 - . . L) are stored, e.g. in memory 42 of the device. 
5 In order to derive accurate and meaningful 

parameters over a wide range of measurement conditions, 
the blood glucose level g as well as the environment tem- 
perature Te is varied during the calibration phase. For 
example, the environment temperature is varied over at 
10 least 5 °C, preferably at least 10 °C, e.g. by carrying 
out indoors and outdoors measurements, and the glucose 
level is varied by at least 100 mg/dl, e.g. by the pa- 
tient having a snack and by delaying and/ or reducing in- 
sulin - 

15 The calibration phase can e.g. extend over 

two days and include at least 10 reference values per 
day. Several reference values should be recorded in the 
periods during which the glucose level and/ or teirperature 
are varied as described above in order to obtain a full 

20 record of these events - 

Alternatively to or in addition to an inten- 
sive calibration phase of two days, an extensive calibra- 
tion can be carried out during a period of e.g. 15 days 
that allows the device to adapt" to a given user. During 

25 this extensive calibration phase, reference measurements 
will again be carried out, e.g. invasively, even though 
at less frequent intervals. 

The data recorded during the calibration 
phase can be used for finding appropriate values for at 

30 least part of the parameters a^. For this purpose, the 
values obtained by function F according to eqfuation (1) 
are compared against the reference values g(ti) or 
against values derived therefrom, and those parameters a^ 
are determined for which this comparison gives a closest 

35 match. 

In a most simple approach, the parameters a^ 
can be obtained from a conventional least-squares fitting 
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algorithm. Suitable algorithms are known to a person 
skilled in the art and are e.g. described by Press, Teu- 
kolsky, Vetterling and Plannery in ^^Numerical Recipes in 
C"/ Cambridge University Press, 2^^ edition, 1992, Chapter 
5 15. For evaluating the function P at the times t^ through 
t]^, only the input values si at the times closest to t^ 
through t^ are required. 

This simple approach, however, will only ex- 
ploit part of the available information- In particular, 

10 it ignores the infoirmation obtained by the measurements 
of the input values Si(t'j) at times t'j other than the 
times ti through t]^. 

In an advanced approach., the reference values 
g(t^) are used to calculate a prediction (interpolation) 

15 of the actual glucose levels at times between the meas- 
urement times ti, t]^, in particular at all times t'l 
... t'L. Then, the deviation of this prediction from the 
value of function F for the corresponding input values s^ 
is calculated and the total deviation is minimized by 

20 varying the parameters a^. 

. An empirical, semi -empirical or theoretical 
model of the variation of the glucose level in a body can 
be used for calculating the prediction (interpolation) . 

An advantageous model is based on the under- 

25 standing that the rate of change of the glucose level is 
limited. For hioman beings, a typical maximum rate of in- 
crease is ffincr - 3.5 mg- dl"l -min""^ and a typical maximum 
rate of decrease is g'decr " ^ mg-dl""^ -min^^ as well. This 
allows to predict a set S of possible glucose values for 

30 any time t between the timesD through t^ as depicted in 
Fig. 4. S is the set of values delimited by lines of 
slope g±ncr 9^ deer extending from the measured points 

g{ti) . 

Taking this model into account, a possible 
35 calibration procedure is based on the following steps: 

Step 1: The patient undergoes the calibration phase 
as mentioned above in which the K reference values 
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g(ti) and the LXN input values sj{t'.j^) are measured 
and recorded. 

- Step 2: Equation (1) is fitted to the measured ref- 
erence values g(t2.) ... g(tK) evaluating 

5 

fi = F(si(ti) ... SN(ti), ao - . . Sij^) (3) 

at each time t^ and comparing f^ to g(ti) . If the 
input values sj (t^) at time t± are not known (be- 

10 cause none of the t'j^ matches t^ exactly, an esti- 

mate of the values sj (t^) from measured input values 
Sj (t'k) for at least one t'k close to t^ can be 
used. Then the parameters ai ... a]y[ are varied to 
find a set of parameters where the total deviation 

15 between the values f^ and g(ti) is at a minimum, 

e.g. by minimizing the sum of the squares of all f i - 
This basic fitting process provides a set of start- 
ing values for the parameters a^ in the following 
step 3 . 

20 - Step 3: The deviation of 



F(t'i) = F(si(t'i) ... SN(t'i), ao ... aM) (4) 

for all times t'^ from the prediction S at the cor- 
25 responding times t'^ is minimized by varying the pa- 

rameters a^. This can e.g. be achieved by defining, 
for each time t'i, a predicted distribution S(t'i) 
of the glucose value and by calculating a deviation 
di by comparing the predicted distribution S(t'i) 
30 with the value F(t'i). In the model of Fig. 4, a 

suitable deviation d^ can e.g. be defined as 



di = 



>(t'i ) - 'Sinax(fi ) if Pit'j_ ) > S^(t\ ) 
S^inax(fi ) - F(t'i ) if Fit'i ) < S^^{t'i_ ) (5) 
0 otherwise 
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with Sinin(t'i) Si^ax^^'i) being the range of the 

set S of Fig. 4 at time t'i, i.e. 

Sinax(t'i) = min(g(tj)+ ff±ncr ' i^' ± ^ t j ) , 
5 fl^decr -(tj+i - t'i)) (6a) 

and 



10 



15 



Smin(t'i) = inax(g(tj)- ffdecx-*(t'i - t j ) , 

g(tj+i)- STincr '(tj+l - t'i)) (6b) 

where t j is the closest of the times ti ... t^ prior 
to t'i. 

The parameters a^ ... ajj can then e.g. be foimd by 
minimizing the value 

D = (7) 



numerically. Corresponding techniques are known to 
the person skilled in the art and e.g. described in 
20 Chapter 10 of the book '^Numerical Recipes in C" 

cited above . 

It must be noted that step 2 is optional if 
the starting values of step 3 are obtained by some dif- 

25 ferent method, e.g. from typical values, or if step 3 
uses an algorithm that does not require starting values 
for the parameters. Alternatively, step 3 can be omitted 
if the results from step 2 are to be used directly. 

It must further be noted that equations (5) 

30 through (7) are advantageous examples but can be replaced 
by other suited definitions. 

For example, instead of using a prediction S 
that gives a simple range, a prediction providing a prob- 
ability density S (g, t'i) can be used, indicating the 

35 probability to observe a given glucose value g at time 
t'i. Such a probability can e.g. be derived from an em- 
pirical or semi -empirical model that predicts how prob- 
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able a given value of the glucose level is at time t'i, 
given the reference values g(tj). Apart from the refer- 
ence values, a suitable model can e.g. take the physio- 
logical parameters of the patient (e.g. body weight) as 
5 well as events during the calibration phase (e.g. food 
intake, insulin administration etc.) into account for im- 
proving the accuracy of the prediction. 

Equation (7) can also be replaced by any 
other suitable measure for the deviation of the fiinction 

10 F from the prediction S. In particular if the probability 
of a certain deviation d^ is known, the formula for D 
should be defined in such a manner that its minimum coin- 
cides with the set of parameters having the highest sta- 
tistical probability. For details, we refer to the book 

15 ''Numerical Recipes in cited above. 

Calibration is preferably carried out with a 
system as shown in Fig. 3, where an external computer 102 
can be. connected to the device 100 through interface 41. 
Computer 102 can instruct device 100 to start a calibra- 

20 tion process, whereupon device 100 can be disconnected 
from the computer and be applied to the patient for car- 
rying out above step 1. Then the reference values g(ti) 
are entered into computer 102, and the measured input 
values Sj(t'i) are transferred to computer 102 via* inter- 

25 face 41 . Above steps 2 and 3 are carried out in computer 
102 and the resulting parameters a^ are transferred back 
to device 100, which, after a final test of the perform- 
ance of the calculated parameters ai, is then ready for 
regular operation. 

50 Even though the capabilities of computer 102 

may be integrated directly into device 100, it is gener- 
ally advantageous to use a separate computer system for 
the convenience of its use and its computational power. 



35 b) Shift correction during- calibration: 

During the above basic calibration, movements 
of the patient or other events may cause device 100 to 
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change its position in respect to the patient's body. 
Displacements of this type will usually lead to a change 
in signal that should be accounted for. 

For taking such shifts into accoiint, it is 
5 advantageous to introduce additional auxiliary parameters 
^00' ^01' • • • ^OP during the above calibration steps. As- 
suming that aQ is a purely additive parameter in function 
F (such as in the example of equation (2)), equations (3) 
and (4) above are replaced by 

10 

fi = F(si(ti) ... S]vj(ti), 0, ai ... a]v[) + 

aoo*^o(ti) + ... + aop-bp(ti) (3') 

and 

15 

F(t'i) = F(si(t'i) ... SN(t'i), 0, ai ... ajyc) + 

aoo-^^O(t'i) + + aop-bp(t'i) , (4') 



where the functions bj_(t) are 0 "unless the time t is in 

20 the range ... ^^+1/ where they are 1. 

In other words, the additive parameter ag of. 
function F is set to 0 (or, equivalently, another fixed 
value) , and it is replaced by parameter aoo in time in- 
terval Tq • • • '^l' parameter agi in time interval Ti . . . 

25 T2' ®tc. 

The times Tq and Xp are the start and end 
times of the calibration phase and the other times are 
the times when a '"shift" of device 100 is detected during 
the calibration phase. Such a shift can e.g. be detected 

30 because at least oxic \jZ LIit= xiij-jLii; values s^ (such as the 
amplitude Aq or frequency fg) changes by more than a 
given threshold value As^ during two consecutive measure- 
ments. Details on how to detect such '"shifts" are dis- 
cussed in the section ""shift correction during measure- 

35 ments" below. 

By using equations (3') and (4') instead of - 
(3) and (4), the parameters slqq ... agp and a^ ... aj^ can 
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be determined using steps 2 and 3 described in the previ- 
ous section. The parameters ai ... a^ can then be used 
during normal operation of the device. 

As to additive parameter slq, that parameter 
5 can be roughly approximated to be the median or average 
of parameters ago • • • ^^OP/ it is preferably deter- 

mined from later recalibration measurements as described 
in section ^Recalibration" below. 

Instead of using additive parameters aoo - • • 
10 aop/ multiplicative parameters might be used for this 

kind of correction as well. In that case, equations (3') 
and ( 4 ' ) should be changed accordingly . 

In more general terms, a compensation of 
"'shifts'' or displacements of device 100 during calibra- 
15 tion can be achieved by replacing at least one of the pa- 
rameters, e.g. aQ/ in equations (3) and (4) by 
p 

i=0 

with bi(t) being 1 for < t < T^+i and 0 otherwise. The 
parameters to be replaced in this way are those parame- 

20 ters that are most sensitive to shifts of the device - 

In most cases, it will be sufficient to apply 
this technique to. the one additive or one multiplicative 
parameter in F. (Definition: A parameter a is additive if 
function f (a, . . - ) can be re-written as a + f ' ( . . - ) with 

25 f ' being independent of a; a parameter a is multiplica- 
tive if function f (a, . . . ) can be re-written as a-f " ( . . . ) 
with f " being independent of a) . 

Normal operation: 

30 After calibration of the device, all or at 

least most of the parameters ag . - - a]y[ are known. In a 
very simple device, such as described in WO 02/069791, 
all parameters can be determined completely during cali- 
bration and then equation (1) can be used for determining 

35 the glucose level from the measured input values s^Ct) in 
regular operation. 
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In the following/ however, some additional 
steps are described that allow to improve the accuracy of 
the device - 

5 a) Recalibration 

After the calibration steps described above, 
all parameters a^ are known if it is assumed that no 
shift correction is necessary, i.e. if it is assimied that 
the device is being held at a fixed position on the pa- 

10 tient's body. 

If a shift of the device against the body is 
to be compensated for, at least one parameter, such as 
the additive or multiplicative parameter a.Q, can only be 
determined inaccurately during calibration because the 

15 device may have been displaced during calibration or be- 
tween calibration and regular measurement. In that case 
it is advantageous to carry out recalibration measure- 
ments during regular operation, e.g. once a day after af- 
fixing the device to the body. 

20 A recalibration measurement consists, in a 

simple embodiment, of a single measurement of the glucose 
level g(to) by conventional means. This glucose level is 
then entered into device 100 with a command to carry out 
recalibration . 

25 When ordered to carry out a recalibration, 

microprocessor 38 finds the solution or optimum agreement 
of 



30 



g(to) = F(si(to) ... si^(to), ao, ai, ... ajs^) (9) 



by varying one of the parameters, usually the additive or 
multiplicative parameter aQ- The parameter found in this 
way is then used for following measurements. 

For solving equation (9), the input values 
35 si(to) ... S]s^(to) may be derived from a single measure- 
ment at time tg or from an average, median or interpola- 
tion value of several measurements around time tQ. Assum- 
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ing that parameters ai to a]^ are known, parameter elq can 
then be calculated e.g. numerically by a root finding al- 
gorithm as known to the person skilled in the art. 

A corresponding recalibration means can e.g. 
5 be iirrplemented as a f iannware program for microprocessor 
38. 



b) Shift correction during normal operation 
As mentioned in the section shift correction 
10 during calibration" above, a movement or shift" of the 
device 100 in respect to the body may cause a chsuage in 
measured signals • Even if all parameters are known from 
calibration or recalibration measurements as described 
above, such a shift may invalidate subsequent measure- 
15 ments . 

To avoid this, microprocessor 38 of device 
100 is advantageously programmed to detect such a shift. 
For this purpose, at least one signal value v(t) can be 
monitored, wherein the signal value v(t) is any value 
20 that is derived directly or indirectly from at least one 
of the input values Si(t) and that shows a characteristic 
shift when device 100 is moved in respect to the pa- 
tient's body. 

In particular, the signal value v(t) ceui be 
25 one of the following: 

One of the input values si(t); for exaitrple, fre- 
quency f 0 or amplitude Aq can be used since both 
these values show a change when device 100 is 
moved . 

30 - The glucose value g derived from function F in 

equation (1) . This value also shows a change when 
the device is moved. 

Any intermediate result generated during the 
evaluation of function F that shows a easily de- 
35 tec ted change when device 100 is moved. 

Fig. 5 shows a typical shift of signal value 
v(t) when device 100 is displaced along the patient's 
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body at a time ts . As can be seen, prior and after the 
event, the signal value is fairly continuous (e.g. lin- 
ear) while there is a sudden change between the measure- 
ments before and after time ts. 
5 To detect a shift of this type, the following 

three steps are carried out at a given time t: 

Step 0 : Calculate an extrapolated signal value 
Vgxt^*^) as an extrapolation from a number of pre- 
vious signal values v. Advantageously, Ve^t^^) is 
10 calculated only from signal values v older than t 

- At. At is a window length, which cem e.g. be 5 
minutes if one measurement is carried out each 
minute . 

- Step 1: Determine an actual signal value Vact(t) 
15 from one or more current signal values v. Advan- 

tageously, Vact(^) is calculated from a median or 
average of the signal values within the time win- 
dow t - At and t. 

Step 2: Compare the actual signal value Vact^^.) 

20 to the extrapolated signal value Vext^^) 

sume that a "shift" has occurred if the values 
differ by a given threshold amount- This thresh- 
old should be larger than a typical noise-induced 
variation between consecutive signals and is e.g. 

25 in the order of 5% of a typical value of v(t) if 

v(t) = fo(t) is used. If the change exceeds the 
threshold amount, a shift correction procedure is 
started- 

The shift correction procedure includes the following 
30 steps: 

Step 3: Define the exact time tg of the shift. 
This can e.g. be done by iterating over a given 
number of recent signal values v(t) , e.g. the 
values in the above time window t - At and t, and 
35 looking for the largest change of consecutive 

values v{tj[) and v(ti_i) . 
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Step 4: Derive a shift correction Av from an ex- 
trapolation of older values (e.g. the extrapola- 
tion Vext(t^) inentioned above) and from values 
measured after the time tg of the shift. For ex- 
5 ample, the difference or ratio between 

- the median or average of the signal values in 
intearval tg ... t emd 

- the extrapolation v^xt^^) 

can be calculated and be used as shift correction 
10 Av. If the difference is used, the shift correc- 

tion will be an additive correction to be added 
to V, otherwise it will be multiplicative correc- 
tion to be multiplied to value v. 
Step 5 : Use the shift correction for correcting 
15 subsequently measured glucose - values . The spe- 

cific implementation of this step depends on the 
definition of the signal value v{t) . Examples: 

If signal value v(t) is equal to an input value 
Si(t), such as fo(t) or Ao(t), subsequently 
20 measured input values should be corrected by 

Si(t) + Av (additive correction) or by Si(t)-Av 
(multiplicative correction) before inserting 
them into frmction F for evaluation. 

- If signal value v(t) is equal to the glucose 

25 value g(t) evaluated from function F, Av can be 

added to or imiltiplied with the returned func- 
tion value. Alternatively, if F has an additive 
or multiplicative parameter aQ, that parameter 
can be corrected by addition of or multiplica- 
30 tion with Av. 

Corrections for other types of signal values v(t) 
can also be implemented by correcting the input 
values, parameters, intermediate results or re- 
turn value of function F in order to make sure 
35 that function F returns the same results before 

and after time tg. 
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Above steps 0 to 5 can be iitplemented in a 
shift correction by suitable firmware in microprocessor 
38 of device 100. In general, the shift correction should 
be able to 

5 - detect a displacement of device 100 along the 

body of the patient, e.g. based on steps 0 to 2 
above or any other method that is able to deter- 
mine a sudden shift in a signal value, 
- determine an effect of the shift on the measured 
10 glucose level, e.g. based on step 4 above, and 

correct the measured glucose levels after the 
shift to compensate for the determined effect. 
It must be noted that the signal value used 
in steps 0 to 3 does not need to be the same as the one 
15 used in steps 4 and 5. It may be advantageous to use a 
raw. input signal, such as fg and Aq for sensitively de- 
tecting a displacement of the device in steps 0 to 3, 
while it may be easier to carry out the correction on the 
function's F return value or an additive or multiplica- 
20 tive parameter ag in steps 4 and 5 . ■ 

Range monitoring: 

As mentioned above, an important purpose of 
device 100 is to provide a prediction of the time when a 
25 patient's glucose level may cross given safety limits. 

For this purpose, microprocessor 38 comprises 
a software- iittplemented predictor that tries to predict 
when, at an earliest time, the glucose level g is likely 
to fall below a lower limit gjoin and/or to rise above an 
30 upper limit or gmax* '^YPicax vaxues for gjixLn the 
order of 50 to 80 mg/dl, e.g. 70 mg/dl, and for gmax they 
are above 160 mg/dl, e.g. 250 mg/dl. 

Such predictors have been known to rely on 
the maximum rate of decrease is g'decr inentioned above, 
35 asstaming that the first derivative g of the glucose 

level will never fall below the maximum rate of decrease. 



wo 2005/053526 



23 



PCT/IB2003/005704 



i.e. gr > -gdecr Q^decr is defined to be a positive 

value) . 

It has been found, however, that this type of 
prediction can be improved. It has been found that not 
5 only the first derivative gr of the glucose level is lim- 
ited, but .also the second derivative gr . Typical lower 
and upper limits gr" and gr**" were found to be both at 0.1 
mg-dl~l •min'"^ . 

This is illustrated in Fig. 6 showing a se- 

10 ries of glucose level measurements g{t) indicated by 

dots. The lines pi and p2 represent worst-case decay pre- 
dictions starting from time tO . pi is calculated on the 
mere assumption that g > -cfdecr ' while p2 is calculated 
from the refined assumption that gr > —ffdecr 3' — * 

15 As can be seen, the time tl where prediction pi reaches 
9min is smaller than the time t2 where prediction pi 
reaches gmin- Hence, using prediction p2 allows to avoid 
unnecessary alerts and allows a more precise prediction. 

To make a prediction of type p2, it is neces- 

20 sary to use not only the actual values g(tO) of the glu- 
cose, but also a first derivative gr (tO) thereof. In the 
example of Fig. 6, time t2 - tO can e.g. be calculated 
from g(tO)/ g (tO) , gr" and gr^ecr rising simple analysis. 

Instead of calculating a time t2. where a 

25 worst-case prediction g(t) is expected to reach gmin' 
is also possible to make a worst-case prediction at a 
time t + At, where At is a fixed ^'safety margin" of e.g. 
20 minutes, and to compare this worst-case prediction g{t 
+ At) e.g. to the lower threshold value gxnin- 

30 worst-case prediction is below the threshold value, an 
alert is issued. 

In general, range monitoring will therefore 
advantageously calculate a prediction of the glucose 
level from an estimate of the current value of the glu- 

35 cose level g(tO) as well as its derivative gr (tO) , taking 
into account that the prediction must fulfil the condi- 
tions gr > -ffdecr ^ ^ "Sr" and/ or gr < cr±^cr ^ ff^ 
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This type of monitoring can be used in the 
device 100 but also in any other type of device that has 
a detector for repetitively measuring the glucose level 
of a living body. The prediction can, in particular, be 
5 used to provide an alert if the worst-case time until a 
hypoglycemia (g(t) < gndn) hyperglycemia (g(t) > gmax) 
is below a given threshold time. 



10 Remarks : 

As it will be clear to the person skilled in 
the art, the methods described above can also be carried 
out with devices different from the one of Figs. 1 and 2, 
such as any of the devices shown in WO 02/069791 (taking 

15 into account that the temperature compensation described 
above will require the addition of .a second temperature 
sensor) . 

Most aspects of the present invention, such 
as a temperature compensation, shift correction and vari- 

20 ous calibration methods, also work with devices using 

other types of sensors, such as optical sensors or induc- 
tive sensors. 

While there are shown and described presently 
preferred embodiments of the invention, it is to be dis- 

25 tinctly xinder stood that the invention is not limited 

thereto but may be otherwise variously eiribodied and prac- 
ticed within the scope of the following claims . 



